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Abstract 

We first review the usefulness of the Wigner distribution functions (WDF), 
associated with Lindblad and pre- master equations, for analyzing a host of 
problems in Quantum Optics where dissipation plays a major role, an arena 
where weak coupling and long-time approximations are valid. However, we 
also show their limitations for the discussion of decoherence, which is gener- 
ally a short-time phenomenon with decay rates typically much smaller than 
typical dissipative decay rates. We discuss two approaches to the problem 
both of which use a quantum Langevin equation (QLE) as a starting-point: 
(a) use of a reduced WDF but in the context of an exact master equation (b) 
use of a WDF for the complete system corresponding to entanglement at all 
times. 
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I. INTRODUCTION 



In recent years, there has been widespread interest in dissipative problems arising in a 
variety of areas in physics. Historically, particularly in quantum optics and NMR (where 
weak coupling and long-time approximations are valid), master equations for the time evolu- 
ation of the density matrix have been extensively utilized in analyzing such problems [1-6] 
and it has also being found useful to re-express such equations in terms of the Wigner dis- 
tribution function [1,4]. However, with the recent burgeoning interest in decoherence, which 
is generally, a short-time phenomenon with decay rates much smaller than typical dissi- 
pative decay rates, a new element arises. Recognizing that conventional master equation 
approaches are often not adequate, some investigators have used path integral methods [3]. 
However, we have found that the simplest and most physically appealing approach to the 
problem is via use of generalized quantum Langevin equations [7], supplemented by use of 
Wigner distribution functions [8]. Our focus is on quantum Brownian motion: a quantum 
particle interacting with a heat bath and moving in an arbitrary external potential. For the 
investigation of decoherence phenomena, we have found [9,10] that it is important to distin- 
guish between two different physical scenarios viz. (a) complete entanglement between the 
quantum particle and the heat bath at all times and (b) the system in a state in which the 
oscillator is not coupled to the bath at, say, t = and such that the bath is in equilibrium 
at temperature T. Thus, it takes a characteristic time of the order of 7 _1 (where 7 is a 
typical dissipative decay time) for the complete coupling to occur and for the whole system 
to come into thermal equilibrium. 

Scenario (b), with the assumption of an uncoupled initial state is that used in the deriva- 
tion of master equations, including the exact master equation of Hu, Paz and Zhang [10]. 
However, this scenario is inappropriate for the calculation of very short decoherence decay 
times and, in general, one must abandon master equation methods for this purpose. For 
analysis of scenario (a), by its nature, it is necessary to abandon master equation methods. 
Instead, we have employed [9] a more general method due to Ford and Lewis [11]. Finally, 
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the presence of a classical external force fit) brings a new dimension to the problem. 

The familiar master equation of quantum optics is in Lindblad form [2], which guarantees 
that the density matrix is always positive definite during time evolution. In the derivation of 
this equation [1,6], rapidly oscillating terms are omitted by the method of coarse-graining in 
time; the high frequencies correspond to the oscillator frequency c^o and, in the usual weak 
coupling limit, c<Jo >> 7, where 7 is a typical decay constant. 

We have referred to the equations obtained prior to coarse-graining in time as pre-master 
(or pre-Lindblad) equations [6,12] and such equations have been used extensively in other 
areas of physics [3,4]; other authors have simply referred to them as master equations but, 
to avoid confusion, we reserve the latter term for equations in Lindblad form. Pre-master 
equations, like the master equation, describe an approach to the equilibrium state. This 
equilibrium state is the same in either case [12], but with pre-master (non-Lindblad) equa- 
tions the approach can be through non-physical states of negative probability. However, 
as we demonstrate here, pre-master equations have other advantages vis-a-vis master equa- 
tions: (a) they lead to the exact expression for the mean value of x(t) (as obtained from the 
exact Langevin equation for the problem), (b) they lead, in the classical limit (h — > 0), to the 
familiar Fokker-Planck equation of classical probability, (c) the exact master equation [23] is 
for long times of pre-master form, positivity being preserved through the time dependence 
of the coefficients. However, the pre-master equations can exhibit different forms, depending 
on whether momentum or coordinate coupling forms of H have been used. 

It is thus clear that there are many facets to the overall question of analyzing dissipative 
systems. Our purpose is to delineate the strengths and weaknesses of the various approaches. 
As will become apparent, one common thread which will permeate the discussion is the 
usefulness of the Wigner distribution for all scenarios. 

As already indicated, there is an important dividing line between long-time (t » 
and short-time (t « 7~ x ) phenomena. Thus, in Sec. II, we discuss the quantum Langevin 
equation for scenario (a), which describes quantum particle and bath entanglement at all 
times, and we consider in turn both the stationery process (corresponding to the absence of 
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an external force fit) in which case we are dealing with the system in thermal equilibrium) 
and then the case of f(t) ^ 0. This naturally leads us to a comparison of the Langevin 
equations of motion with those derived (in the Ohmic limit and for an oscillator) using 
the conventional master equation techniques of Quantum Optics and NMR. Thus, as a 
starting-point, we write down in Sec. Ill, the quantum Langevin equation for the Ohmic 
model. In addition, we present, in the interaction representation, the two forms (momentum 
coupling and coordinate coupling) of the weak-coupling pre-master equations as well as the 
corresponding master equation (which is the same for both couplings [12]). 

Then in Sec. IV, we consider in detail the momentum and coordinate coupling forms, 
respectively, of the pre-master equations as well as the master equation. In Sec. V, we discuss 
the corresponding Wigner quantum probability distributions and their classical limit. 

To incorporate short-time phenomena, a more sophisticated analysis is required, requiring 
both the Langevin equation discussed in Sec. II as well as a generalization to a Langevin 
equation for the initial value problem, which is the theme of Sec. VI. First, we present the 
solution of this equation in terms of the solution of the stationery equation, the corresponding 
Green function and the fluctuation force. Next, we write this solution in a form local in 
time. Moreover, we point out that the equivalent master equation, as distinct from those 
discussed in Sees. Ill to V, is exact in that the coefficients are time-dependent for which 
explicit expressions are readily determined using the aforementioned solution of the Langevin 
equation. Then, in Sec. VII, with a view toward calculation of quantum probability functions 
for superposition states and the corresponding decay of the interference term (decoherence) 
and developing a framework which will accomodate many different possibilities, we are 
led again to a consideration of Wigner distribution functions, both for the whole system 
of quantum particle and bath and also (after averaging over the bath variables) for the 
reduced system. We also show that it is sometimes useful to express the solution for the 
Wigner function at time t in the form of an integral of a transition probability acting on the 
initial Wigner function; in particular, we are led to a generalization of the classical Kramers 
equation, in the sense that it is quantum, applies to an arbitrary heat bath and allows for 
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a time-dependence in the external force. In Sections VI and VII, we confined ourselves to 
f(t) = and, next, in Sec. VIII, we show how these various results can be applied to the 
calculation of decoherence, both for the case of (a) entanglement at all times and (b) the 
initial value problem where the system and bath are initially decoupled. In addition to our 
detailed discussion of decoherence in coordinate space, we also comment on results obtained 
for decoherence in both phase space and momentum space. The case of fit) ^ is also 
considered. 

Figs. 1 and 2 present an overview of the paper which should hopefully enable the reader 
to "see the wood from the trees." In particular, in Fig. 1, we see that specification of the 
complete Hamiltonian H leads immediately to a quantum Langevin equation, the exact form 
for which depends on whether one is considering entanglement at all times or an initial value 
problem. From this equation, one is led to Wigner quantum probability distributions from 
which one may obtain, in particular, probabilities in either coordinate space or momentum 
space. These are the tools necessary for consideration of various problems such as decoher- 
ence phenomena. Fig. 2 deals with the weak coupling limit of exact master equations and 
delineates their shortcomings. 

II. GENERALIZED QUANTUM LANGEVIN EQUATION: ENTANGLEMENT 

FOR ALL TIME 

In recent years, there has been widespread interest in dissipative problems arising in 
a variety of areas in physics. As it turns out, solutions of many of these problems are 
encompassed by a generalization of Langevin's equation to encompass quantum, memory, 
and non-Markovian effects, as well as arbitrary temperature and the presence of an external 
potential V(x). As in Ref. [7], we refer to this as the generalized quantum Langevin equation 
(QLE) 
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f dt'n(t - t')x(t r ) + V'(x) = F{t) + f(t), 



(2.1) 



where V'(x) = dV(x)/dx is the negative of the time-independent external force and fi(t) is 
the so-called memory function. F(t) is the random(fluctuation or noise) force and f(t) is a c- 
number external force. In addition, it should be strongly emphasized that "- the description 
is more general than the language -" [7] in that x(t) can be a generalized displacement oper- 
ator (such as the phase difference of the superconducting wave function across a Josephson 
junction). 

A detailed discussion of (2.1) appears in Ref. [7]. In particular, it was pointed out 
the QLE corresponds to a macroscopic description of a quantum system interacting with 
a quantum-mechanical heat-bath and that this description can be precisely formulated, 
using such general principles as causality and the second law of thermodynamics. We also 
stressed that this is a model-independent description. However, the most general QLE can 
be realized with a simple and convenient model, viz., the independent-oscillator (10) model. 
The Hamiltonian of the 10 system is 



Here m is the mass of the quantum particle while rrij and Uj refer to the mass and frequency 
of heat-bath oscillator j. In addition, x and p are the coordinate and momentum operators 
for the quantum particle and qj and pj are the corresponding quantities for the heat-bath 
oscillators. 

The infinity of choices for the rrij and Uj give this model its great generality [7]. In 
particular, it can describe nonrelativistic quantum electrodynamics [7], the Schwabl-Thirring 
(S-T) model [13], the Ford-Kac-Mazur (FKM) model [14], and the Lamb model [15]. 

In this context, it should be noted that, whereas H in (2.2) has been put into a form 
in which all the heat-bath oscillators interact with the central oscillator of interest but not 



H 




(2.2) 
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with each other, we have shown that this is the most general H one can write down to 
describe most types of dissipation encountered in the literature; in particular, it is unitarily 
equivalent to the FKM model in which all the oscillators are coupled (as discussed in Sec. 
VE of Ref. [7]). Also, of interest, in connection with our later discussions, is the velocity- 
coupling model in which the coupling is through the particle momentum (as distinct from 
the particle coordinate as in (2.2)), corresponding to the Hamiltonian 

Hvc = ^ (p + E + V(x) + £ + \m^) q ]^ - xf(t). (2.3) 

As shown explicitly in [7], H vc is obtained from H by means of two unitary transformations 
which change the p, pj, and qj but which do not change the coordinate x. We have refrained 
from putting primes on the changed variables since the context will always ensure that no 
confusion will arise. However, for future purposes, we note that, in proceeding from (2.2) to 
(2.3), 

P^P-J2 m j UJ j ( lj- ( 2 - 4 ) 

j 

Thus, it is not surprising that both the 10 model and the vector coupling model lead to the 
same identical GLE (2.1). Morover, since both Hamiltonians (2.2) and (2.3) are unitarily 
equivalent to the Hamiltonian of non-relativistic quantum electrodynamics (QED) we see 
the close analogy with the equivalence of the x ■ E and p ■ A couplings in QED. Whereas 
most of our previous work has used the 10 form of the Hamiltonian, we found it convenient 
in [6] to use the VC form for the derivation of a non-Markovian master equation. 

Use of the Heisenberg equations of motion leads to the GLE(2.1) describing the time 
development of the particle motion, where 



ji{t) = rrijU? cos(u jt)6(t), 



(2.5) 



is the memory function, with 6{t) the Heaviside step function. Also 

F(t) = J2m J ^(t), (2.6) 

3 

is a fluctuating operator force with mean (F(t)) = 0, where q h {t) denotes the general solution 
of the homogeneous equation for the heat-bath oscillators (corresponding to no interaction). 
These results were used to obtain the results for the (symmetric) autocorrelation and com- 
mutator of F(t), viz., 

\(F(t)F(t') + F(t')F(t)) 
1 f°° 

= - / duoRe\jl(uj + i0 + )]frw coth(hw/2kT) cosb(t - if)] (2.7) 
n Jo 

Oft roc 

\F(t), F(t')} = — / dwRe{jl(u + i0 + )}u 

VK JO 

Xfflnw(t-t'). (2.8) 
Here fl(z) is the Fourier transform of the memory function: 

jliz) = / dtfi(t)e izt . (2.9) 
Jo 

Equation (2.7) is referred to by Kubo [16] as the second fluctuation-dissipation theorem 
and we note that it can be written down explicitly once the GLE is obtained. Also, its 
evaluation requires only knowledge of Rep,(u). On the other hand, the first fluctuation- 
dissipation theorem is an equation involving the autocorrelation of x(t) and its explicit 
evaluation requires a knowledge of the generalized susceptibility a(cu) (to be defined below) 
which is equivalent to knowing the solution to the GLE and also requires knowledge of both 
Re/2(cj) and lmjl(ou). This solution is readily obtained when V(x) = 0, corresponding to the 
original Brownian motion problem [17]. As shown by FLO [7], a solution is also possible 
in the case of an oscillator. Taking V(x) = \Kx 2 = \muolx 2 , these authors obtained the 
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solution of the Langevin equation (2.1) in the form 

x(t)= f dt'G(t-t'){F(t')+f(t>)} : 
where G(t), the Green function, is given by 



(2.10) 



1 r°° 

G(t) = — dcua(u + iO + )e- tuJ \ (2.11) 

27T J-oo 

with a(z) the familiar response function 



—mz l — izn(z) + K 
It is often convenient to write (2.10) in the form 



x(t) = x s (t) + x d , (2.13) 

where Xd is the "driven" contribution [18] due to the external force fit) and x s is the contri- 
bution due to the fluctuation force F(t). Here we have introduced a subscript s to emphasize 
that x s (t) is a stationary operator-process, in the sense that correlations, probability distri- 
butions, etc. for this dynamical variable are invariant under time-translation (t — > t + to). 
In particular, the correlation, 

C Q (t - t') = ^{x s (t)x s (t') + x s )x s (t)) 

H f°° JlLU 

= -/ du;lm{a{uj + i0 + )} coth — — cosu(t-t'), (2.14) 
7r Jo 2kT 

is a function only of the time-difference t — t'. Furthermore, [18] 

C d (t,t') = ^{x(t)x(t f )+x(t')x(t)) 

= C (t-t') + (x(t))(x(t')), (2.15) 

where (x(t)) is the steady mean of the driven motion. 



Also, taking the Fourier transform of (2.10), we obtain 



x(u)=a(u){F(u) + f(u)}, (2.16) 

where the superposed tilde is used to denote the Fourier transform. Thus, x(u>) is the Fourier 
transform of the operator x(t): 



x{uj) = / dtx{t)e iu}t . (2.17) 
We have now all the tools we need to calculate observable quantities. 

III. QUANTUM LANGEVIN EQUATION FOR THE OHMIC MODEL AND 
CORRESPONDING WEAK COUPLING MASTER EQUATIONS 

As a special case of (2.1), we see that in the particular case of an oscillator with spring 
constant K, the Langevin equation for Ohmic coupling (constant 7), has the form 



mi + m-fx + Kx = F(t) + f(t), (3.1) 

Moreover, 

<*(*)> = <p(f)> = - 7 <P(*)> " K(x(t)) + /(f). (3.2) 

m 

These are exact results. In the weak coupling limit, the corresponding pre-master and 
master equations for the density matrix may be obtained. First, a word about notation [6]. In 
the Schrodinger representation, the density matrix for the whole system (quantum oscillator 
plus heat bath) is denoted by p whereas p denotes the reduced density matrix, obtained 
after a partial trace has been carried out over the heat bath variables (which, in effect, 
eliminates the very rapid stochastic fluctuations). The equation for p is what we refer to as 
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a pre-master equation and it includes high frequency oscillating terms exp{±2icu t} , where 
u is the oscillator frequency. If the equation for p is time-averaged (which is equivalent to 
discarding the high-frequency oscillating terms) then the resulting equation is called a master 
equation for the time-averaged density matrix which is denoted by p. The corresponding 
quantities in the interaction representation are a, a and a. Explicitly [6], in the interaction 
representation, 



a = cxp 



{fM-f}, (3.3) 



a(t)=Tr R {a(t)}, (3.4) 
(where Ttr refers to the trace over the reservoir), 

— = — + exp{±2iu t} terms, (3.5) 

and we have corresponding equations in the Schrodinger representation. 

Moreover, the form of any of these density matrix equations depends on whether one is 
using the momentum or the coordinate coupling forms. In the interaction representation 
[6,12] the momentum coupling form can be written 



^ = -^{^ooth^[p(t)Ap(t),a]]- i [p(t),x(t)a + ax(t)]}, (3-6) 
while the coordinate coupling form can be written 



it = ~2^ {m ^ oCOth 2^ [;r(t) ' M^+^WMt^ + vpit)]}- (3-7) 
From either of these equations one obtains the same master equation by discarding 
the terms explicitly oscillating at frequency 2uo (often referred to as the rotating-wave 
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approximation (RWA)) and this is best achieved by transforming from the variables x,p to 
a and a + [6]. 

The corresponding master equation is [12] 



and it is a simple matter to verify that this has the Lindblad form of the master equation 
familiar in quantum optics [1,5,6]. On the other hand, outside this community, the pre- 
master equations have been used extensively by Caldeira and Leggett [3] and others. We 
wish to emphasize that all three of the above equations lead to the same equilibrium state: 
p cq = exp{—Ho/kT}; i.e., detailed balance is obeyed. The difference is in the approach 
to equilibrium. For the pre-Lindblad equations this can be through (unphysical) states in 
which p is not positive definite. Note that this form of the equilibrium state holds only if 
Hq is the free oscillator Hamiltonian (given below in (4.3)). 

As already alluded to, whereas the Lindblad form of the master equation has the merit of 
ensuring a positive definite matrix element at all times (in contrast to pre-master equations) 
it has the disadvantage of being subject to the RWA. On the other hand, there are undesirable 
features associated with the non-positivity associated with the pre-master equations [19] but 
it should be remarked that all master equations give an incorrect form for very short times, 
as can be shown most simply by calculating the mean-square displacement, defined in (8.8) 
below, by means of both an exact calculation and a master equation approach. On the 
other hand, after a short initial time of the order of the "bath relaxation time" (zero for 
the Ohmic case, r for the single relaxation time case) the exact equation has the pre-master 
form. This can be seen from [10], where in equation (5.4) we gave explicit expressions for 
the Ohmic case or in equation (A5) for the single relaxation time case. Another point is 
that the pre-master equation gives the exact equation of motion for the mean motion, the 
master equation does not. 



da 
~dt 



1_ 

Ah 



{ coth 




+i([x,pa + up] - [p,xa + ax})}. 



(3.8) 
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IV. MEAN VALUES FOR THE VARIOUS FORMS OF THE DENSITY MATRIX 
A. The Momentum Coupling Form of the Pre-Master Equation 

In the Schrodinger representation (3.6) becomes [6] 



dp 1 7 r 1 fkvo , . , r _ _ . . 

-rrr = Tz[H d ,p] ~ T^i coth 7^UflP> )P,P]]-Wxp + px\}, 4.1 

at in in muo 2kl 

where Ha, the Hamiltonian for the driven oscillator, is given by 



H d = H -xf(t). (4.2) 

and where 



H.= £ + \K* (4-3) 

is the free oscillator Hamiltonian [18]. 

The corresponding equations for the mean of x and p are 



d (n) d 

-rAx) = — - l(x), -tM = -rru^Kx) + f(t). (4.4) 
at m dt 

If we eliminate (p) between these two equations, we get 



^ + /M + ^ {x) = m (4 . 5) 

This, of course, is just what we would get from the quantum Langevin equation. 
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B. The Coordinate Coupling Form of the Pre-Master Equation 

In the Schrodinger representation (3.7) becomes 



-7^ = jfi[ H d,p\ ~ ^ muJ ° coth ^kf^ x ' +i[x,PP + PP}}- (4.6) 

Here, and also below, we neglect the small energy shifts because they are of no conse- 
quence as far as the present discussion is concerned. The corresponding equations for the 
mean of x and p are 

d (n) d 

Jt (x) = ^, Jt (p) = -rruv 2 (x) - 7 <P> + /(*)■ (4-7) 

If again we eliminate (p) between these two equations, we get the same equation (4.5) 
as in the momentum-coupling case. In other words, both forms of the pre-master equation 
lead to what we would get from the exact quantum Langevin equation. 

C. The Master Equation 

In the Schrodinger representation (3.8) becomes 

dp 1 7 frw 1 

dt = ih [Hi ' P] ~ 4^ {c ° th 2^T Ml +™bfr 

+i([x,pp + pp] - [p, xp + px])}. (4.8) 
The corresponding equations for the mean of x and p are [18] 

^ = ^-i<»>,^ = ^>-^ + /W. (4.9) 
Eliminating (p) between these two equations, we get 

d 2 (x) d(x) . 9 7 2 . . . f(t) 

-i/ + T^ + M + ^)W = ^. (4.10) 
14 



Because of the spurious (7 2 /4) term, this is uot the equation one gets from the quantum 
Langevin equation. In fact, this equation is identical with the equation obtained from use 
of an RWA approximation to the exact Hamiltonian [20]. 



V. QUANTUM PROBABILITY DISTRIBUTIONS 

These are most simply obtained by using the Wigner function which is defined by the 
relation [8] 

1 Z" 00 11 11 

W(q,p; t) = ^ due^/ h p(q -- q +-- t ). (5.1) 

The translation from the operator equation for p to the quasi-classical equation for the 
Wigner function is most easily carried out by use of the Bopp operators [8], as a result of 
which we get the correspondence 

PP^{p- y ^)W, pp^(p+ J-^)W. (5.2) 

It follows that the corresponding equation for the Wigner function is 

dW p dW , dW r . ,dW 
at m dq op op 

+ | {(1 + A) ^ + (1 _ A) ^ } 



, L |\ ^ h d 2 W „ ^mfku d 2 W\ 

+l( * + 2»ft ¥ +(1 " A| Vv I' (5 ' 3) 



where A = 0, 1 and -1 correspond to the master, momentum coupling pre-master and 
coordinate coupling pre-master equations, respectively, and where 



2(JV + h = coth(hcu /2kT). (5.4) 



As a check, we note that the equilibrium distribution, 
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^•') = wW^ { -|£iS; } ' (5 ' 5) 

satisfies this equation for any value of A if f(t) = 0. It is also instructive to write (5.3) 

explicitly for each value of A. 

(a) Momentum Coupling Pre-Master (A = 1) 

dW p dW 2 dW r . .dW 
+ mu q— f(t)- 



dt m dq dp dp 

dqW , Ar L h d 2 W 
dq 2 mujQ dq 



+7^ + # + t^' (5-6) 



(b) Coordinate Coupling Pre-Master (A = —1) 

dW p dW 2 dW „ .dW 

~aT = JT + m ^o9^ /W^ - 

dt m dq dp dp 



dpW , AT 1. d 2 W 
^_ + 7 (JV + -)m^— . (5.7) 



(c) Master (A = 0) 



<9W p 2 dW r/ s dW 
~74T = 7T~ + mw o9^ 

7 <9gW <9pW 
/ 1\ f ft <9 2 W mhujo d 2 W) , , 

We note that, in the high-temperature limit, (5.7) reduces to the well-known classical 
Kramers equation [21], generalized to include time- dependence in the external force. Actu- 
ally, (5.6) is not inconsistent with (5.7) since the values of the momentum p appearing in 
both equations are not the same, as discussed in the Introduction (while emphasizing that 
q is the same in both equations). On the other hand, (5.8) does not reduce to the correct 
classical result. Earlier discussions of different forms of the master and Fokker-Planck equa- 
tions may be found in [4] and [22]. In the next section, we will show how the exact master 
equation is obtained from the solution of the initial value quantum Langevin equation. As 
above, it is found to be desirable to express the result in the form of an equation for the 
corresponding Wigner distribution (see (6.7) below). 
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VI. QUANTUM LANGEVIN EQUATION: INITIAL VALUE PROBLEM 



In this and the following section, we will take f(t) = 0. The Langevin equation for the 
oscillator with given initial values is given by [10] 

mi + /* dl/fiCt - t')x(t') +Kx = -fji(t)x(0) + F(t), (6.1) 

JO 

and the general solution is given by 

x (t) = mG(t)x(0) + mG(t)x(0) + X(t) (6.2) 
where we have introduced the fluctuating position operator, 

X(t) = t dt'G(t-t')F(t'). (6.3) 
Jo 

If we assume that at t — the system is in a state in which the oscillator is not coupled to 
the bath and that the bath is in equilibrium at temperature T, we find that the correlation 
and commutator are the same as those for the stationary equation. Also, the Green function 
is the solution of the homogeneous equation, 



mG + f dt'/iH - t')G(t') + KG = 0, (6.4) 
Jo 

with the initial conditions 



G(0) = 0, G(0) = — . (6.5) 

777- 

Typically, the memory function /i(t) falls to zero in a very short time r, called the 
relaxation time of the bath. For times long compared with this bath relaxation time, the 
extra term on the right hand side of (6.1) vanishes, but only for much longer times, times long 
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compared with the oscillator decay time, will this equation become the stationary equation, 
with the lower limit on the integration taken to be — oo. It was then possible to show [10] 
that one could write the Langevin equation (6.1) in the form of an equation that is local in 
time with time-dependent coefficients: 



x + 2T( y t)x + Q 2 ( y t)x = — F(t), (6.6) 

m 

where explicit expressions for Y(t) and Q(t) were obtained in terms of G(t). Furthermore, 
it was shown that these results constitute in essence a derivation of the HPZ exact master 
equation [23] with explicit expressions for the time- dependent coefficients. 

The HPZ equation is perhaps best displayed in the form of its corresponding Wigner 
function equation. Using (5.2), we obtained [10] 

dW 1 dW „ 2/ , dW _ .dpW 

-~r = — p^t + rnn 2 (t) q — + 2r t -g— 

at m oq Op op 

d 2 W d 2 W 

+ hmr(t)h(t)^ + hr(t)f(t)—, (6.7) 

where f2 2 (t), 2T(t), h(t), and f(t) are time-dependent parameters for which we obtained 
explicit expressions in terms of X(t), F(t) and G(t). Since these are rather cumbersome, 
we simply refer to equations (2.19) and (3.8) of [10]. It should also be noted that equations 
equivalent to the HPZ equation have been obtained in [24-26] and in [27,28] for the particular 
case of an Ohmic bath. 

The solution of this exact master equation will be addressed in the next section. However, 
it is prudent to point out a serious caveat at this stage: a serious divergence arising from the 
assumption of an initially uncoupled state is found to be due to the zero-point oscillations 
of the bath and not removed in a cutoff model. As a consequence, worthwhile results for the 
equation can only be obtained in the high temperature limit, where zero-point oscillations 
are neglected [10]. 



18 



VII. WIGNER DISTRIBUTIONS FOR THE INITIAL VALUE PROBLEM 



Our explicit solution (6.2) of the Langevin equation, together with use of Wigner func- 
tions enables us to construct the general solution of the exact master equation. Since the 
latter is an equation for the reduced density matrix, it was natural to proceed by writing 

W(q,p;t) = J dq_J dpW systcm (q,p;q,p;t). (7.1) 

Here W systcm is the Wigner function for the system of oscillator and bath, with q = (qi, q 2 • • ■ ) 
and p = (pi,P2 • • • ) the bath coordinates and momenta. Also, since the initial state for the 
HPZ equation is a product state, it follows that the corresponding Wigner function is of the 
form 



W S ystem(?,P;q,P;0) = W{q, p\ 0) Y[ Wj (q 3 ;,Pj) . (7.2) 

i 

Here, on the right W(q,p; 0) is the initial Wigner function for the oscillator and the product 
is the Wigner function for the bath, in which Wj(qj,pj) is the Wigner function for a single 
oscillator of mass rrij and frequency u>j, 

1 



w Mil Pi) = 



irhcoth(foujj/2kT) 

y ( Pi+™phi ) {73) 

\ rrijhuJjCoth.(hujj/2kT) J 



After some algebra, we were able to write the solution in the form of a transition operator 
acting on the initial Wigner function, 

/oo /*oo 
dq' / dp'P(q,p;q',p';t)W(q\p';0). (7.4) 
-oo J ~ oo 

where P(q,p; q',p'; t), the transition probability, can be written 



19 



P(q,p; q',p';t) 



27rv / detA 

where we have used a dyadic notation with 

/ 

A(t) = 



exp]--R-A- 1 .R[, 



(7.5) 



m 2 X 2 



f(XX + XX 



v f (XX + XX 



(X 2 



R(t) 



(7.6) 



v- {P^)) 

9 -<?(*)> 

Here, in R, the quantities (q(t)) and (p(t)) correspond to the mean of the initial value 
solution with initial values q' and p'. We refer to [10] for further details. We also note that 
transition probabilities are also discussed in [29]. 

In fact, the expression (7.4) for the transition probability is formally the same as that 
for the classical Kramers equation [21]. The difference is that the Green function and 
the mean square of the fluctuating position and velocity operators here are for a quantum 
oscillator interacting with an arbitrary heat bath, while in the classical solution of the 
Kramers equation they are for a classical oscillator interacting with an Ohmic bath. Also, 
the classical solution of the Kramers equation corresponds to taking the stationary limit. 

Finally, it should be emphasized that the key results we obtained in [10], some of which 
will be discussed below, depended only on our solution of the initial value Langevin equation, 
in the course of which we showed equivalence with the HPZ equation. 



VIII. DECOHERENCE 

The problem of decoherence in quantum systems has been of considerable recent interest 
[30-35,9,10]. Decoherence refers to the destruction of a quantum interference pattern and 
is relevant to the many experiments that depend on achieving and maintaining entangled 
states. Examples of such efforts are in the areas of quantum teleportation [36], quantum 
information and computation [37,38], entangled states [39], Schrodinger cats [40], and the 
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quantum-classical interface [41]. For an overview of many of the interesting experiments 
involving decoherence, we refer to Refs. [39] and [42]. 

Much of the discussion of decoherence has been in terms of a particle moving in one 
dimension that is placed in initial superposition state (a Schrodinger "cat" state) corre- 
sponding to two widely separated wave packets (see (8.3) below) separated by a distance 
d. Decoherence is said to occur when the long-time interference pattern is destroyed. The 
key questions asked are, first, under what conditions does decoherence occur and, second, 
what is the decoherence time. In broad outline, the problem is that of a quantum particle 
in a superposition state and interacting with the environment and two distinct physical sce- 
narios present themselves depending on whether the particle is initially decoupled from the 
environment or, alternatively, in thermal equilibrium with the environment at the time it is 
put into the initial state by a measurement. We will now examine each situation in turn, 
using the results developed in Sections VI and VII. 

A. Initial Decoupling of Particle and Environment 

This has been the situation which has examined by most previous investigators and the 
tools of choice are either the Feynman- Vernon influence functional techniques [30] or master 
equation techniques [31-33,10]. It has now become clear that, in general, the Lindblad 
and other weak-coupling master equations are not adequate because they cannot properly 
handle the short times associated with decoherence phenomena [10]. We have recently given 
an exact treatment of the problem based on our exact solution of the Langevin equation for 
the initial value problem and concomitantly, the exact solution of the HPZ master equation 
[10]. As discussed in Sec. VII, the solution is best presented as a Wigner function for the 
reduced density matrix. We recall that the HPZ equation is a master equation with time- 
dependent coefficients for a harmonic oscillator interacting with a linear passive heat bath of 
oscillators. The equation is exact and general within the assumption that in the initial state 
the bath is in equilibrium and not coupled to the oscillator. This assumption of a decoupled 
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initial state is common to all derivations of a master equation, going back at least to the 
work of Wangness and Bloch [43] , who phrased it as an assumption that at any instant of 
time the system is approximately decoupled. Indeed, such an assumption is essential for 
the introduction of the notion of partial trace, i.e., the trace over states of the uncoupled 
bath, key to the existence of any master equation. The bath at t = is in equilibrium at 
temperature T whereas, for the oscillator, we have the freedom to decide its initial state. 
More often, the initial temperature of the bath is taken to be zero or T but it is important 
to emphasize that, depending on the choice, different predictions ensue [10]. Hence, we 
will examine both situations in turn. In both cases, we will need as a starting-point the 
probability density at time t, given by 



where W(x,p, t) is the reduced Wigner function. This enabled us to obtain the general 
result [10] 



This is the only case where one can start with a pure Schrodinger superposition ("cat") 
state, undisturbed by temperature or environmental effects (albeit difficult to achieve in 
practise). When this state at temperature zero is suddenly coupled to the heat bath at 
temperature T, the initial time dependence is then dominated by the "warming up" of 
the particle, which occurs on a time scale of order the decay time 7 -1 . However, since 
7 _1 is generally very much greater than the decoherence decay time r^, decoherence will 






(8.2) 



where W is the Fourier transform of the Wigner function. 



(i) Particle initially at zero temperature 



22 



have occured before the particle has reached the bath temperature, as a result of which an 
important contribution to r d is missed. To see this explicitly, we considered an initial state 
corresponding to two separated Gaussian wave packets. The corresponding wave function 
has the form 



1>(x, 0) = p (exp{-^— |^ } + exp{-^), (8.3) 

where d is the separation and a is the width of each packet. Next, using (8.1), we obtain 
an expression for the probability distribution which consists of three terms, the first two of 
which correspond to the individual wave packets whereas the third is an interference term, 
which consists of a cosine term multiplied by a time-dependent factor. The attentuation 
factor a(t) is the ratio of the coefficient of the cosine term divided by twice the geometric 
mean of the first two terms. We find, for an Ohmic bath, high T and t « 7 _1 [10] 



where ( = rwy. If we suppose that a is negligibly small, we find a(t) = exp{— t/r d } where 
Td = ^!f d 2 ■ This, except for a factor of 6 is exactly the decoherence time that often appears 
in the literature [32,44] for the off-diagonal terms of the density matrix (as distinct from 
the spatial probability discussed here, which is of more interest because its a measurable 
quantity). On the other hand, if the a term in the denominator dominates, then decoherence 
occurs on a time scale ~ t 3 . But, as we have seen above, this result corresponds to a particle 
in an initial state that is effectively at temperature zero, which is suddenly coupled to a heat 
bath at high temperature. The result is therefore unphysical in the sense that the initial 
state does not correspond to that envisioned when we speak of a system at temperature T. 
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(ii) Particle initially at the bath temperature T 



In this case, when the particle and bath are suddenly coupled, the whole system is im- 
mediately in equilibrium at temperature T. In order to analyze this scenario, we simply 
generalize to a state corresponding to a particle with a random velocity v due to the tem- 
perature environment. Next, we obtain the corresponding Wigner function which we average 
over a thermal environment. This leads to the result [9,10] 



— t 2 

a T (t)^exp{ s d 2 }, t«7~\ (8.5) 

W 1 8(a* + a 2 ^t 2 + £?t 2 ) V } 

which for very short times is of the form a-r(t) — exp{— t 2 /rj}, where the decoherence time 
is 

r d = — , (8.6) 

in which v = \J kT/m is the mean thermal velocity. In fact, as we have shown [34,35], the 
latter result can be obtained solely within the framework of elementary quantum mechanics 
and equilibrium statistical mechanics. In particular, since is independent of the param- 
eters of the heat bath, the result given in (8.5) exhibits "decoherence without dissipation" 
[9,34,35]. Thus, very different results are obtained than for case (i). Moreover, it is clear 
that a similar procedure can be carried out for the case of a particle initially in a very general 
state different from that of zero temperature. 

However, as mentioned previously, all of these calculations are subject to a serious caveat 
connected with the assumption of an initially uncoupled state viz. we found [10] that a 
serious persistent divergence related to the zero-point oscillations of the bath and which are 
not removed in a cutoff model. As a consequence, worthwhile results for the equation can 
only be obtained in the high temperature limit, where zero-point oscillations are neglected. 
Such problems do not arise in case (B), which we will now examine. 



24 



B. Particle and Environment Entangled at all times 



In order to describe a state of the system that is entangled at all times, including the 
initial time, it is necessary to abandon master equation methods. Some time ago, a more 
general method applicable to such systems was described by Ford and Lewis [11]. In their 
method, a system in equilibrium is put into an initial state (e.g., a wave-packet state) by a 
measurement and then at a later time is sampled by a second measurement. This method of 
successive measurements has recently been applied to obtain exact results for the problems 
of wave packet spreading and decoherence [9]. For the decoherence problem, one obtains 
the result [9] 

^ = ^{-^k}- < 87 > 

for a free particle initially in the state given by (8.3), where s(t) is the mean square dis- 
placement. For f(t) — 0, we have 

s(t) = {[x s (t) - x s (0)f 

7T 

is the mean square displacement for the stationary process and where the variance is given 
by 



COO 

\ ^lm{a(cu + ?0 + )}coth— — (1 -cosut), (8.8) 
jo 2kT 



wHt) = a > + s(() _ M!W0)]!. (8 . 9) 

The calculations leading to these results were based on the quantum probability func- 
tions, introduced by Ford and Lewis [11], and which may be shown to be closely related 
to Wigner distribution functions. In the high temperature limit, (8.7) reduces to the same 
result given in (8.5). This is not very surprising since (8.5) corresponds to the particle and 
bath both being at the same temperature, albeit uncoupled, at t — (case A(ii)). How- 
ever, when they are suddenly coupled, they immediately become entangled, as distinct from 
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case A(i), where the particle temperature is initially zero so that it takes a time ~ 7 -1 for 
entanglement to occur. 

In our view, this scenario is the most physically realistic one because at all times the 
particle is coupled to the bath. Moreover, we have obtained exact results, for all t, which 
incorporate both arbitrary temperature and arbitrary dissipation and, in particular, contains 
non-trivial t dependences in both numerator and denominator, especially in the case where 
kT Ki hq. Also, even at T = 0, writing a(t) = exp{— b(t)}, we find that b(t) ~ — t 2 log7r for 
the single relaxation time model with (t/r) « 1. [10,47]. 

While our calculations of a(t) have concentrated on the time decay of the interference 
term in the probability distributions, they were derived from Wigner distributions or closely 
related quantities [9,10] from which other measures of decoherence are readily calculated, 
leading to similar conclusions. 

In order to calculate decoherence in phase space, we start by considering any quantum 
state consisting of two identical components separated by a distance d. Then, it is not 
difficult to show that the corresponding Wigner distribution function, W^(q,p) at time 
t = 0, say, is given by 

W {2) (q, p, 0) = N I W (q + ^ p, oj + W (q - ± p, oj + 2 cos ^ W(q, p, 0) j , 

(8.10) 

where W(q, p, 0) is the Wigner function for one of the pairs at t — and N is a normalization 
factor. Applying this result to the two-Gaussian superposition, we concluded that there is 
no decoherence in either phase space or momentum space i.e. decoherence is only manifest 
in coordinate space [45]. 

The effect of an external classical force f(t) on decoherence has also been considered 
[46]. In the absence of dissipation and for negligibly low temperature (which approximates 
the experimental conditions discussed in [42]) we obtained [46] 
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where 



Sd(t) = (x 2 d (t)) 



Jo Jo 



/ dt' dt"G(t-t')G(t-t")g(t' -t"), 



(8.12) 



is the mean-square displacement due to the driving force with 



0(f-O = </(O/(O>- 



(8.13) 



In the case of a random delta-correlated force so that 



g (t'-t")=g5(t'-t"). 



(8.14) 



where g is time-independent, we obtained, for small times characteristic of decoherence 
phenomena, 



In summary, we have shown that the Wigner distribution function is an invaluable tool 
for analyzing dissipative problems in quantum mechanics including investigations dealing 
with decoherence. Also, with regard to the latter, we demonstrated that very different 
results may ensue depending on the choice of whether the particle is initially decoupled 
from the environment or, alternatively, in thermal equilibrium with the environment at the 
time it is put into the intial state by a measurement. 
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